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ABSTRACT 

IRAC imaging of a 4'7 x 4'7 area in the Lockman Hole detected over 400 
galaxies in the IRAC 3.6//m and 4.5//m bands, 120 in the 5.8yum band, and 80 in 



^ ! the S.O^m band in 30 minutes of observing time. Color-color diagrams suggest 

(^ ! that about half of these galaxies are at redshifts 0.6 < z < 1.3 with about a 

^ ! quarter at higher redshifts {z >1.3). We also detect IRAC counterparts for 6 of 

Tj- ! the 7 SCUBA sources and all 9 XMM sources in this area. The detection of the 



counterparts of the SCUBA sources and galaxies at z > 1.3 demonstrates the 
ability of IRAC to probe the universe at very high redshifts. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 



Oh' 
O 

(^ ! Survey — galaxies: Mid- Infrared 



X. 

H . 1. Introduction 

Much effort has gone into deep optical imaging with hope of detecting galaxies at high 
redshifts (Koo & Kron 1992; Ellis 1997). Obtaining redshifts for faint, optically selected 



^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA02138 

^Spitzer Science Center, Caltech, 1200 E. Cahfornia, Pasadena, CA91125 

■^Department of Astrophysics, Oxford University, Keble Rd, Oxford, 0X1 3RH, U.K. 

''Steward Observatory, University of Arizona, Tucson, AZ85721 

^Institut d'Astrophysique Spatiale, bat 121, Universite Paris Sud, F-91405 Orsay Cedex, France 

^Max-Planck-Institut fiir Astronomic, Konigstuhl 17, D-69117 Heidelberg, Germany 

^Subaru Telescope, National Astronomical Observatory of Japan, 650 North A'ohoku Place, Hilo, HI 96720 



-2- 



galaxies require long observations on >8-nieter class telescopes, but most of the galaxies 
identified are at intermediate redshifts (z~l) (Cowie et al. 1996; Cohen et al. 1999; Cimatti 
et al. 2002). Applying faintness criteria to visible- wavelength (UBVRI) surveys is not an 
efficient way of selecting high-redshift galaxies unless color selection methods such as the 
Lyman break technique (Steidel et al. 2003) or extreme red color selection (Thompson et al. 
1999) are applied. Furthermore, optical studies of high redshift galaxies are heavily affected 
by dust extinction because rest-frame UV photons are being observed (Steidel et al. 2003). 

By contrast, far-infrared and sub millimeter surveys are efficient tools for generating 
samples of high-redshift starburst galaxies. For example, SCUBA on the JCMT has made 
it possible to carry out deep extragalactic surveys in the submm band. Now it is known 
that most SCUBA sources are ultra luminous infrared galaxies (ULIRG) at high redshifts. 
They are detected in the submillimeter band, because the infrared luminosity peak at about 
100-200 /xm shifts to the submillimeter and millimeter band at z>2 (Chapman et al. 2003; 
Ivison et al. 2002; Small et al. 2002). 

Little is known about the properties of mid-infrared selected samples, although in prin- 
ciple they should have many advantages. Mid-infrared observations of high-redshift galaxies 
detect light emitted at 1-2 /xm near the peak of the stellar emission. At these wavelengths 
the emission is only modestly affected by dust extinction. The difficulty is that a space obser- 
vatory is essential for sensitive mid-infrared observations. The Infrared Space Observatory 
(ISO) (Kessler 1996) carried out the first deep infrared extragalactic surveys and detected 
galaxies at redshifts up to z = 1.5 in the mid-infrared bands (Oliver et al. 1997; Elbaz et al. 
1999; Aussel et al. 1999; Genzel & Cesarsky 2000; Dole et al. 2001; Sato et al. 2003). With 
its combination of high sensitivity and high spatial resolution in four mid-infrared bands, the 
Infrared Array Camera (IRAC) (Fazio et al. 2004) on the Spitzer Space Telescope (SST) can 
carry out deeper and larger surveys and should be able to detect galaxies at 2; > 3 (Barmby 
et al. 2004). 

This Letter reports the results of a SST survey carried out during the SST In-Orbit- 
Checkout (IOC) period. Observations were taken as one of the Early Release Observations to 
demonstrate the sensitivity of Spitzer and its capability of detecting galaxies at high redshifts. 
The program included both IRAC and MIPS observations of the same Lockman Hole field. 
Egami et al. (2004) and Le Floc'h et al. (2004) present the results of the MIPS observations 
of the same field. The survey area centered on RA(J2000)=10:51:56, DEC(J2000)=57:25:32 
has been covered by surveys at many other wavelengths, including SCUBA, XMM, and deep 
optical/NIR imaging. Of seven SCUBA sources in the survey area, IRAC detected six in 
the 3.6 and 4.5 yum bands and four in the 5.8 and 8.0 ^m. bands. IRAC also detected mid- 
infrared counterparts of all 9 XMM sources. Analysis of the Spitzer counterparts of the 
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SCUBA, XMM and MAMBO sources are presented by Egami et al. (2004), Alonso-Herrero 
et al. (2004) and Ivison et al. (2004). 



2. Observations and data reduction 

The observations^ covered one field-of-view in all four IRAC bands. The observations 
used the 30-s frame time and 25 dithers per pointing, giving an exposure time (i.e., effective 
on-chip charge accumulation time after accounting for readout overheads) of 670 s in each 
band. We used a shorter frame time and many dithers to be able to reject cosmic-rays 
effectively. Because of the dithering, the area covered to full depth was 4'7 x 4f7. 

Data processing began with the Basic Calibrated Data (BCD) delivered by the Spitzer 
Science Center (SSC). These data include flat-field corrections, dark subtraction, and linear- 
ity and flux calibrations. Additional steps included pointing refinement, distortion correction, 
and mosaicking. Cosmic rays were rejected during the last step by sigma-clipping. 

Source extraction used DAOPHOT because the images are so crowded in the 3.6/xm and 
4.5/im bands even for the short observations reported here. Given the I'.'S FWHM of the 
point spread function (Fazio et al. 2004) in both bands, almost all faint sources are point-like. 
We used a signal/noise ratio of 2.5, and a PSF with full-widths at half-maximum (FHWM) 
of I'.'S, I'.'S, 2'.'0, and 2'.'0 in the 3.6, 4.5, 5.S, and S.O |Um bands respectively as parameters for 
DAOPHOT. 

Star /galaxy separation was based on morphology in deep R-band images taken with the 
Suprime camera on the Subaru Telescope and on a {V — I vs I — K) color-color diagram 
(Huang et al. 1997; Wilson et al. 2001). We identified 13 stars in the field with K < 20.5. 
Though fainter stars may be missed, their density is extremely low (Glazebrook et al. 1994). 
Table 1 gives the numbers of galaxies and stars detected in each band along with the limiting 
flux. The number of galaxies detected in the MIPS 24/im band (Egami et al. 2004; Le Floc'h 
et al. 2004) is also in Table 1. Fewer sources were detected in the two long wavelength bands 
because of the high zodiacal background at those wavelengths. All objects detected at 5.S 
or S.O /xm are detected at both 3.6 and 4.5 yumThere are about 25000 "beams" in the field, 
and the 3.6 /xm source counts thus correspond to 57 beams per source. 



^Spitzer Program ID 1077, executed 2003 November 20 



-4- 



3. The Properties of Galaxies Detected by IRAC 

With no redshifts available for most of the IRAC galaxy sample, we rely on color-color 
diagrams to estimate redshifts. In the rest near-infrared wavelength range, the most impor- 
tant spectral feature is the 1.6 jum bump, which is caused by the H~ opacity minimum in the 
atmospheres of cool stars. This feature is seen in the spectral energy distributions (SEDs) of 
all galaxies except AGN-dominated ones and has long been considered a photometric redshift 
indicator (Wright et al. 1994; Simpson & Eisenhardt 1999; Sawicki 2002). kt z > 0.6, the 
bump shifts beyond the K band, and /j,(3.6/im) > fi^{2.2fim). For z > 1.3, the bump passes 
3.6 /xm, and /j,(4.5yum) > /,^(3.6/im). Local, dusty starburst galaxies such as M82 also have 
red colors with /,^(4.5/im) > f^{3.6fim), and we need to require /,^(3.6yum) > /j,(2.2yum) to 
identify them. 

Figures 1 and 2 show the IRAC data along with the tracks followed by typical galaxy 
SEDs as redshift increases, and Table 2 summarizes the observed galaxy number distribution 
as a function of color, which imply redshift. Figure 1 shows that the observed colors are 
consistent with the model tracks. The most-populated redshift bin in Table 2 for all four 
IRAC bands is 0.6 < z < 1.3. This bin includes about 50% of all detected IRAC galaxies. 
Figure 1 shows a few sources with very red K-[3.6] color. Wilson et al. (2004) found that these 
sources are actually the same extreme red objects as those detected in optical-near-infrared 
surveys using the (R-K)>5 color selection. 

The color-color diagrams will not give accurate redshifts for AGN, which do not exhibit 
the 1.6/im bump in their SEDs. X-ray observation is a good way of detecting AGN. This 
field has 9 XMM sources (Hasinger et al. 2001), whose colors are plotted along with other 
galaxies in Figures 1 and 2. These 9 sources probably represent most of the AGN in the field 
and constitute only a few percent of the galaxies detected. Thus they should not affect our 
statistics on redshift distribution. 

Unlike the XMM sources, whose color distribution fills the entire color space in Figures 1 
and 2, most of the SCUBA source counterparts have colors indicating z > 1.3. Using the 
deep radio counterpart positions to identify the SCUBA sources in the optical and Spitzer 
images, Egami et al. (2004) did a detailed analysis of SED for the SCUBA sources. They 
found that LE 850. 1^, LE 850.4, LE 850.7, and LE 850.14 have the 1.6yum bump in their 
SEDs. This implies that their SEDs are dominated by stellar light, and thus their IRAC 
colors following the galaxy color tracks well in Figures 1 and 2 with expected redshfits z > 2. 
This is consistent with most of the SCUBA sources being dusty galaxies at z > 2 (Chapman 
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et al. 2003; Egami et al. 2004). LE 850.8 is the only object for which we have both SCUBA 
and XMM detection. Egami et ah (2004) found that LE 850.8 has a power-law continuum. 
Its IRAC colors are not close to any of the tracks in Figures 1 and 2. Another SCUBA source, 
LE 850.18, also has colors not close to any of the tracks and a power-law SED (Egami et al. 
2004). In the deep R-band image, however, there are more than 2 objects within 3" radius 
centered on the radio counterpart positions for both LE 850.8 and LE 850.18. IRAC will 
not be able to resolve them if all components emit in the IRAC bands. Hence we suspect 
foreground contamination in the IRAC flux densities in both cases and are cautious in using 
IRAC colors to classify both objects. 

There are 32 sources in this fleld detected in the MIPS 24 ^m band at or above the 
5-sigma confldence level. Their IRAC counterparts show that half of the MIPS sources are 
aX z > 1.3. Given the 5a limiting flux of 200 /iJy, these MIPS sources at 2; > 1.3 have a 
minimun infrared luminosity of 10^^ Lq, thus are LIRGs and ULIRGs (Le Floc'h et al. 2004). 

The extragalactic background light (EBL) from both resolevd and unresolved extra- 
galactic sources is an indicator of the total luminosity of the Unverse (Madau & Pozzetti 
2000). Because the intermediate redshift galaxy population contributes most to the EBL 
(Elbaz et al. 2002), integral resolved galaxy light (IGL) in a deep survey should give a good 
lower limit on the EBL. Table 3 shows the results of integrating the fluxes, uli,, from IRAC 
counterparts of objects identifled by SCUBA (Scott et al. 2001), XMM (Hasinger et al. 
2001), and MIPS 24 /im images (Egami et al. 2004; Le Floc'h et al. 2004). The IGL at the 
5.8 //m and 8.0 yum are 1.31 and 1.02 nW m^^ sr~^, very close to the ISO LW2 (6.75 yum) 
IGL, 1.7±0.5 nW m"^ sr^^ (Ahieri et al. 1999; Elbaz et al. 2002). Given that the XMM 
sources in this fleld are neither stars or clusters, it is safe to assume that most of them are 
AGN (Lehmann et al. 2000; Alonso-Herrero et al. 2004). The light contributed by the XMM 
sources is constant in all four IRAC bands because of the flat SED of AGN in the mid- 
infrared. Because the total light decreases with wavelength, the XMM source contribution 
increases from 3.7% at 3.6 fim. to 13.6% at 8 /im. This is consistent with the 15±5% XMM 
source contribution found by the ISO survey (Fadda et al. 2002). As shown in Figure 2, most 
of the MIPS 24 yum sources are at intermediate and high redshifts. Most of them must be 
LIRG's or even ULIRG's (Egami et al. 2004). The starburst light percentage increases from 
20% at 3.6 /xm to 62% at 8 yum. AGN and starburst galaxies together contribute about two 
thirds of the IGL at 8 /xm. The bulk of the light at 3.6 yum and 4.5 /xm appears to come from 
evolved stars in galaxies, while the 8 yum and 24 yum light traces the star forming regions in 
galaxies. This explains the large contribution of the 24yum selected sample to the 8//m IGL. 
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4. CONCLUSIONS 

Even with only 30 minutes observing time, IRAC detected 419 galaxies in a 4'.7 x 4'7 
area in the Lockman Hole. Based on their colors, most of the detected galaxies are likely 
to be at intermediate redshifts 0.6 < z < 1.3, but >25% appear to be at z > 1.3. 30% of 
the 5.8 /im sources and 50% of the 8.0 /im sources show colors indicating higher redshifts 
(z>1.3). The counterparts of the SCUBA sources have colors consistent with galaxies at 
z > 1.3. The integrated galaxy light seen in our sample is consistent with 1.7±0.5 nW m^^ 
sr^^ derived using the ISO survey data at the 6.75 /im. The AGN fraction, as revealed by 
XMM, increases from 3.7% in the 3.6 /im band to 13.6% in the 8.0 /im band. The MIPS 
population, most of which are starburst galaxies at intermediate redshifts, contributes up to 
60% in the 8.0/im band. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through Contract Number 
1256790 issued by JPL. 
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Fig. 1. — Color-color diagram for objects detected at 3.6 yum and 4.5 /im. Filled circles are 
objects also detected at 24 /im, open circles are remaining sample. Arrows indicate objects 
not detected at K, points are plotted at 3a" limits. Blue diamonds denote XMM sources, 
and red circles denote SCUBA sources. Solid lines show color tracks for elliptical galaxies 
(blue), Sbc galaxies (green), and M82 (red). Expected colors for ellipticals and spirals are 
based on rest frame SED's from Coleman et al. (1980) in visible light and Lu et al. (2003) 
in the infrared. M82 tracks are derived from the model of Granato et al. (2000) and ISO 
observations (Sturm et al. 2000). The tracks are derived by shifting the SED without adding 
any evolution; numbers marked along each track indicate redshift. In this diagram, galaxies 
in the upper left quadrant are likely to be at 0.6 < z < 1.3, those in the upper right are 
likely to be at z > 1.3, and galaxies in the two lower quadrants are likely to have z < 0.6 
with starburst galaxies to the right. 
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Fig. 2. — Color-color diagram for galaxies detected at 8 /xm. Blue diamonds show XMM 
sources, red dots SCUBA sources with numbers according to Scott et al. (2001), and black 
dots are objects detected at 24 /xm. Open circles show the remaining sample galaxies. Color 
tracks are the same as those in Figure 2. 
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Table 1 
Number of Objects detected in the IRAC bands 



Band 


Galaxies 


Stars 


Total 


Mag 


Limit(5o-,AB) 


3.6Mm 


419 


13 


432 




23.73 


4.5/im 


403 


13 


413 




23.77 


5.8/im 


120 


5 


125 




21.90 


8.0^111 


80 


5 


85 




21.68 


24.0^m 


32 





32 




18.15 



Table 2 
Galaxies in each IRAC Color Bin 



Kai 



[3.6], 



[3.6], 



[4.5]yis Estimated Redshift 3.6/im 4.5^m 5.8/im 8.0^in 24/im 



> 

> 

< 

< 






z > 1.3 


113 


122 


35 


24 


16 





0.6 < z < 1.3 


191 


176 


57 


34 


8 





z < 0.6 


101 


81 


27 


22 


8 





z < 0.6 


14 


24 


1 









Table 3 
Integrated Galaxy Light u I,, (nW m^^ sr^^) in the IRAC Bands 





3.6^m 


4.5/im 


5.8/^m 


8.0^m 


All galaxies detected 


3.47 


2.43 


1.31 


1.03 


SCUBA counterparts 


0.05 


0.05 


0.05 


0.04 


XMM counterparts 


0.13 


0.12 


0.13 


0.14 


MIPS counterparts 


0.68 


0.52 


0.43 


0.62 



